nx=300;
ny=300;
width=5;
distance=30;
start_45=2;
start_m45=17;


leftw=10;
rightw=10;
bottomw=50;

rands=[5,8,3,12];
[xx,rand_num]=size(rands);



%========UNIFORM WIDTH=============================
% geo1=ones(nx,ny);
% for xnum=-nx+start_45:distance:nx
% for i=1:nx
%     for w=1:width
%         x=i;y=i+xnum+w;
%         if ((x>0) && (x<=nx) && (y>0) && (y<=ny))
%             geo1(x,y)=0;
%         end
%     end
% end
% end
% 
% for xnum=start_m45:distance:2*nx
% for i=1:nx
%     for w=1:width
%         x=i;y=-i+xnum+w;
%         if ((x>0) && (x<=nx) && (y>0) && (y<=ny))
%             geo1(x,y)=0;
%         end
%     end
% end
% end
%=================================================





%================RANDOM WIDTH=======================
geo1=ones(nx,ny);
for xnum=-nx+start_45:distance:nx
    width=rands(ceil(rand*rand_num));
for i=1:nx
    for w=1:width
        x=i;y=i+xnum+w;
        if ((x>0) && (x<=nx) && (y>0) && (y<=ny))
            geo1(x,y)=0;
        end
    end
end
end

for xnum=start_m45:distance:2*nx
width=rands(ceil(rand*rand_num));
for i=1:nx
    for w=1:width
        x=i;y=-i+xnum+w;
        if ((x>0) && (x<=nx) && (y>0) && (y<=ny))
            geo1(x,y)=0;
        end
    end
end
end
%====================================================

%==============ADD BOUNDARIES===============
[lx,ly]=size(geo1);
bufdat=zeros(leftw,ly);
geo1=[bufdat;geo1];
bufdat=zeros(rightw,ly);
geo1=[geo1;bufdat];
[lx,ly]=size(geo1);
bufdat=zeros(lx,bottomw);
geo1=[bufdat,geo1];
%===========================================




